Stochastic equation for the erosion of inclined topography 
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We present a stochastic equation to model the erosion of topography with fixed inclination. The 
inclination causes the erosion to be anisotropic. A zero-order consequence of the anisotropy is the 
dependence of the prefactor of the surface height-height correlations on direction. The lowest higher- 
order contribution from the anisotropy is studied by applying the dynamic renormalization group. 
In this case, assuming an inhomogenous distribution of soil material, we find a one-loop estimate of 
the roughness exponents. The predicted exponents are in good agreement with new measurements 
made from seafloor topography. 
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The rich complexity of the Earth's surface, both on 
land and beneath the sea, is the result of physical mech- 
anisms ranging from tectonic motion to surficial erosion 
jj],!]. Despite this variation, however, geologic surfaces 
show a certain degree of universality: they may often be 
characterized as self-affine |^^] over some range of length 
scales. This means that, if h(S,to) is the height of the 
surface at position x at some time to, then the "rough- 
ness", measured by the height-height static correlation 
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grows as x 



function C(x) = ((h(x,t ) - h(Q,t )) 2 } 
where a is called the roughness exponent [Q. Empirical 
measurements of a are numerous. While many indicate 
that a is small (0.30 < a < 0.55) H|), a number of other 
measurements show it to be large (0.70 < a < 0.85) 
Moreover, some measurements indicate that a 
crosses over from large to small values as length scales 
become greater than approximately 1 km Q . Motivated 
by these findings, we propose that the large values of a at 
small length scales may be explained by the influence of a 
preferred direction — downhill — for the flux of eroded ma- 
terial. We derive an anisotropic noisy diffusion equation 
to describe erosion at the small length scales where the 
preferred direction is fixed throughout space. Under the 
additional assumptions that the flux of eroded material 
increases with increasing distance downslope and that the 
dominant effects of noise arc fixed in space, we find, using 
the dynamic renormalization group (DRG), a first-order 
estimate of the roughness exponents. New measurements 
of our own, made from the topography of the continen- 
tal slope off the coast of Oregon, are in good agreement 
with our predictions. We find that our anisotropic the- 
ory significantly enriches previous isotropic continuum 
models ||,[lCj for two reasons. First, it predicts that cor- 
relations differ in different directions, and and second, 
it predicts that these correlations decay quantitatively 
differently than they do for isotropic topography. 

Figure [l] depicts the framework for our theory: a sur- 
face h on a two-dimensional substrate. We refer to h 
generically by the term landscape, and note that it in- 
clination is fixed. The unit vector is the "growth" 




FIG. 1. Schematic configuration of an anisotropic land- 
scape for the case d = 2. 

direction, which is measured downwards from the top of 
the slope. The preferred, downhill, direction is given by 
the unit vector e||, while ej_ represents a vector perpen- 
dicular to e || and e^. Later, when applying the DRG, 
we will generalize to landscapes on a d-dimensional sub- 
strate; in this case e_L represents the subspace of all di- 
rections perpendicular to eii and e^, and has dimension 
d — 1. The configuration is completed by selecting fixed 
boundary conditions at the top of the slope, x\\ =0, or 
by imposing the symmetry x\\ — > —x\\. 

Due to the preferred direction eii in Fig. [|, the sta- 
tistical properties of h may be anisotropic. Thus, if h 
is self-affine, we expect different roughness exponents for 
correlations measured in each of the directions eu and 
ej_. Thus we define at\\ and aj_ such that C||(iE||) ~ a;?" 
for correlations along a fixed transect = const., and 
C±(x±) ~ x°^- for correlations along a fixed transect 
x?i = const., where in general a\\ ^ a±. These relations 
can be summarized in the single scaling form 



C(x l{ , x ± ) - 6 Q » C^-^y , b~^x±), 



(1) 



where Q is the anisotropy exponent. The exponents ay 
and aj_ are related through a±_ — a||/£||. The exponent 
C|i accounts for the different rescaling factors along the 
two main directions. Since the space is anisotropic, when 
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performing a scale change, we must rescale xu and x± 
by different factors 6|| and b±, respectively, if we are to 
recover a surface with the same statistical properties. We 
assume in our model that £|| = log fej^ / log 6|| = const. 

We seek a single stochastic equation for the land- 
scape height h. Whereas others || have advocated 
the now classical, isotropic, non-conservative interface 
growth equation due to Kardar, Parisi, and Zhang (KPZ) 
jnj, we assume here that the underlying soil is locally 
conserved such that 

d t h = -V-J + r), (2) 

where J is the current of soil per unit length. The soil 
however is not globally conserved, since it is lost at the 
bottom boundary. We also allow local conservation to 
be broken by the addition of a stochastic noise term 77, 
discussed below. 

Physically, the current J is expected to reflect two ef- 
fects. First, we expect a local isotropic diffusing com- 
ponent, tending to smooth out the surface. Second, we 
expect an average global flow of dragged soil, directed 
mainly downhill. Thus we postulate the following form 
for the current: 

J= -uVh-Wwh. (3) 

The first term corresponds to Fick's law for diffusion, 
and represents the isotropic relaxational dynamics of the 
soil. The second term represents the average flow of soil 
that is dragged downhill, either due to the flow of wa- 
ter or the scouring of the surface by the flow of the soil 
itself. The direction of this term is given by the vector 
V||ft = 9|ihen. The term T plays the role of an anomalous 
anisotropic diffusivity. In order to gain insight into the 
role of r, consider the case in which erosion results from 
the stress exerted on the soil bed by an overland flow q 
of water, where q is the volumetric flow rate though unit 
area perpendicular to the direction of steepest descent. 
The greater q is, the stronger is the stress More- 
over, since q flows downhill, it increases with distance 
downslope. Thus T must be an increasing function of 
xii. Since the fixed inclination implies that h increases 
with 111, we choose to parameterize the anomalous dif- 
fusion as a function of the height such that T = T(h) 
||. Defining T(h) = A + g(h), with g(0) = and 
G(h) = f g(h)dh, we substitute Eq. (||) into Since 
g(h)d\\h — [dG(h)/dh]d\\h — d\\G(h), where we have used 
the chain rule for the second equality, we obtain 

d t h = v^dpi + is ± V 2 ± h + d^Giti) + rj, (4) 

where v±_ = v and v\\ = v + Ao. 

Even in the absence on any nonlinearity, fundamen- 
tal conclusions may be drawn from (H) . By setting g = 
(i.e., by considering T(h) = Ao = const.), we obtain a lin- 
ear equation which is an anisotropic counterpart of the 



Edwards- Wilkinson equation M. In can then be easily 
shown [jf) that the correlation functions along the main 
directions ey and are inversely proportional to the 
square root of the diffusivities m and v±_ respectively, 
that is, Cx/Cii ~ (Vn/^jJ 1 / 2 . In other words, since 
the preferred direction gives > the topography 
is quantitatively rougher, at all scales and by the same 
factor, in the perpendicular direction than in the parallel 
direction. 

In order to obtain more information on the scaling 
properties of Eq. (|]), we have studied it using the DRG. 
Assuming that T(h) is an analytical function, we can per- 
form a Taylor expansion in powers of h. Since all odd 
powers of h must vanish in order to the preserve the joint 
symmetry h — > —h, J — > —J in Eq. (||), we are left at 
lowest order with g(h) ~ A2/1 2 . By dimensional analysis 
one can check that all the terms in this expansion are 
relevant under rescaling. However, the flux Q(x\\) of the 
erosive agent (water or soil) flowing on the surface should 
grow no faster than Q(x\\) ~ x d . Then, taking h ~ x\\, 

we find that the terms in g(h) should be of order h d or 
less. Specializing to the case of d = 2 (i.e., real surfaces), 
we then find it reasonable to truncate g at second order, 
such that Eq. (||) takes the form 

d t h=v ll djh + v x Vlh+^df ] (h s ) + r l , (5) 

where A = A2. Note that Eq. (||) differs from the 
anisotropic driven diffusion equation of Hwa and Kar- 
dar |f|| because the form of our current J is suggested 
not only by symmetry arguments, but also by the physics 
of erosion. 

We now address the issue of noise. We distinguish 
two different sources. First, we may allow a term of 
"annealed" noise, rjt(x,t), depending on time and po- 
sition, and describing a random, external forcing, due to, 
for example, inhomogeneous rainfall. We assume that 
this noise is isotropic, Gaussian distributed, with zero 
mean, and uncorrelated such that (rjt(x,t)r]t(x' ,£')) = 
2D t 8^ d \x — x')S(t — t'). Second, we may have a term 
of "quenched" noise to account for the heterogeneity of 
the soil, mimicking the variations in the erodibility of 
the landscape j8|. We represent this randomness by a 
source of Gaussian static noise r) s {x), with correlations 
(rj s (x)r] s (x')) = 2D s 5^ d \x — x'). This form of noise has 
been previously proposed to model soil heterogeneity in 
cellular automata models of fluvial networks JlJ] . In the 
following we consider the limits (i) r/ s — (D s <C Dt), 
corresponding to a situation of random external forcing 
and homogeneous composition of soil, and (ii) rj t = 
(D s ^ D t ), representing the limit in which the external 
forcing is constant and the most essential source of noise 
is the inhomogeneous composition of the soil. 

Application of the DRG follows the procedure used in 
Refs. (T^JT^]. In Fourier space we proceed by integrating 
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over the shell of large wave vectors Ae~' < k < A, where 
A is the wave vector upper cutoff and e l is the rescaling 
factor, and by subsequently rescaling the system back to 
its original size through the transformation x± — > e l x±, 
xh — > e'^ x xii, h — > e la± h, and i — > e iZi t. The anisotropy 
is explicitly included in the exponent £j_ = Ci7 • To 
lowest order in perturbation theory, both limits (i) and 
(ii) above provide the same form for the renormalization 
group flow equations: 

— = \{z± + 2a±_ - 2Ci_ - 

^1 = Di(KiZ± - 2aj_ - Cl — d + 1), 
dl 

where i = i, s stands for the limits (i) and (ii) 
above, respectively. Here X; is an effective cou- 
pling constant, depending on the type of noise: 
A* = AAifd-iA d - 2 /2^ /2 ^ /2 in (i), and A s = 

\D s K d ^ d - A l2v^ 2 v Z l 2 in (ii), with K d = S d /(2Tr) d 
and Sd the surface area of a d-dimensional unit sphere. 
The value of the correction factor m is Kt = 1 and k s — 2. 
The flow equations for vy and Di are exact to all orders 
in the perturbation expansio n Jl3| , ^6[ . They provide us 
with the exact result zj_ = 2 Jl7[. The effective coupling 
flows under rescaling as 

^ = Xi(si - 3Ai), (6) 

where = d c — d, and d c is the critical dimension 
for each particular limit, d c — 2 and d c = 4. The 
stable fixed points of (||) are A* = for d > di l) and 
A* = Si /3 for d < d c . For d > d c ^ the critical exponents 
attain in both limits their mean- field values a± = 0, 
Cl IF = 1, and zf F = 2. On the other hand, for d < 
the critical exponents computed at first order in the e- 
expansion are: 

a± W = §' C±(i)=l + ^r- (7) 

The physically relevant dimension for erosion is d = 2. 
In the limit of thermal noise this corresponds to the 
critical dimension. By continuity, the exponents are 
a_L = aii = and £j_ = C|| = !• This result is consis- 
tent with a flat landscape, with logarithmic corrections 
to the roughness 0. However, we still expect anisotropy 
to appear in the prefactor of the correlation functions Cm 
and Cj_, as argued above. On the other hand, in the 
limit of static noise we are below the critical dimension, 
and (0) is applicable. Substituting e s = 2 we obtain the 
roughness exponents 

a± = \ c~ 0.83, ay = £i = ~ ~ 0.63. (8) 



The values (g) predicted for a± and an are in reason- 
able agreement with previous measures made at small 
length scales However, these measurements were 

either averaged over all directions or the direction of the 
measurements was not reported. Thus, to check our re- 
sults with a natural landscape that has an unambiguous 
preferred direction, we have analyzed digital bathymetric 
maps of the continental slope off the coast of Oregon. In 
this case the slope results from the relatively abrupt in- 
crease in the depth of the seafloor as the continental shelf 
gives way to the deeper continental rise. Figure 2a shows 
one portion of this region. Here the main feature of the 
topography is a deep incision called a submarine canyon. 
In this region, submarine canyons are thought to have re- 
sulted from seepage-induced slope failure which oc- 
curs when excess pore pressure within the material over- 
comes the gravitational and friction forces on the sur- 
face of the material, causing the slope to become unsta- 
ble. Slope instabilities then create submarine avalanches, 
which themselves can erode the slope as they slide down- 
wards. 

Figure 2b shows the height-height correlation functions 
C || and Cj_, corresponding, respectively, to the parallel 
and perpendicular directions of the seafloor topography 
in Fig. 2a. The computation of C± follows from its defini- 
tion but the computation of Cn requires some comment. 
The fluctuations measured by C|| must be defined with 
respect to an appropriate average profile. Briefly, one 
expects that geologic processes other than erosion (e.g., 
tectonic stresses) are responsible for long-wavelength de- 
formation in the parallel direction. We may estimate such 
systematic corrections by computing the mean profile in 
the parallel direction: h av (xn) = LJ 1 J dx±h(xu, x±). 
We then compute Cm from the fluctuations of the de- 
trended surface h — h— h av (x\\). ^From both Cy and C± 
we find that the least-squares estimates of the roughness 
exponents, a« ~ 0.67 and a± ~ 0.78, exhibit a surpris- 
ingly good fit to our theoretical predictions (|j). 

We have also measured Cy and C± in some desert en- 
vironments. In these cases (not shown), we did not ob- 
tain conclusive power law scaling, but we always found 
C±/C\\ > 1, as predicted by the linear theory. Thus, 
while the example of Figure 2 may be in some sense spe- 
cialized, one of our main predictions — that the topog- 
raphy in the perpendicular direction is rougher than the 
topography in the parallel direction — seems to be of fairly 
general validity. 

In conclusion, we note that the main elements of our 
theory are the conservation of the eroded material, ran- 
domness of either the landscape or the forcing, and the 
presence of a preferred direction for the material trans- 
port. The latter assumption leads to an anisotropic equa- 
tion that applies, in principle, to any erosive process with 
the appropriate lack of symmetry. In the usual geological 
setting, however, the anisotropy applies specifically to a 
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FIG. 2. (a) Digital map of a submarine canyon off the 
coast of Oregon, located at coordinates 44°40' N, 125°45' 
W. The vertical axis represents depth below sea level. Dis- 
tances are measured in meters, (b) Height-height correlation 
functions computed along the parallel (Cm) and perpendicu- 
lar (C±) directions. Solid lines are least-squares fits to the 
scaling region. 



surface of fixed inclination which, in turn, implies that 
our theory should only apply locally, to the relatively 
small scales where the preferred direction of transport is 
approximately constant. Because the anisotropy should 
vanish at large length scales, these large scale features 
should be presumably described by an isotropic theory, 
such as the KPZ equation P,p3|. Indeed, the KPZ equa- 
tion predicts exponents that are approximately consis- 
tent with large scale observations. Since these predictions 
differ from ours, it may be possible to use our results to 
distinguish statistically between features of the landscape 
due to erosion and features due to larger-scale processes, 
such as tectonic deformation. 
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